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Abstract 



We propose a scheme to measure quantum Stokes parameters, their fluctu- 
ations and correlations. The proposal involves measurements of intensities 
and intensity-intensity correlations for suitably defined modes, which can be 
produced by a combination of half wave and quarter wave plates. 
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Polarisation of light is a well studied subject. Text book (l| descriptions of light polarisa- 
tion are based on the Stokes parameters. Consider, for example, a light beam propagating in 
z direction and decomposed into two components polarised along the orthogonal directions 
e,-, i = 1, 2: 



ikz — iwt 



E{z ) t) = {ha l + t 2 a 2 )e ikz - iwt + c.c. (1) 
The four Stokes parameters characterising such a beam are defined by |I[ 

O I 1 2 i I 1 2 o I |2 I |2 

oq = |«i| + \a 2 \ ; oi = |oti| — | ai2| ; 

^2 = (a 1 «2 + «i«2) > S3 = — 2(a*«2 — aic^X (2) 

which, as is well known, lead to the representation of light polarisation by a point on the 
Poincare sphere. All states of linear polarisation lie in the xy plane, whereas north and 
south poles represent states of circular polarisation. All other points on Poincare sphere 
represent states of elliptic polarisation. Methods to measure all four Stokes parameters are 
also well known. 

The Stokes parameters have been generalised in several ways 0-0 . For example Wolf 
has shown how the polarisation characteristics of partially coherent light can be calculated. 
In such a case one has to average quantities defined by (2), over an ensemble distribution 
of ol\ and«2. More recently, there has been considerable interest in generalising these ideas 
to discuss the polarisation characteristics of a non classical field possessing very strong 
quantum fluctuations In such a case the field E in (1) becomes an operator with aii 

and a 2 replaced by the annihilation operators a\ and a 2 . We can then define quantum Stokes 
operators f7|, in analogy to (2), as : 

S = (a\ai + a\a 2 ), (3a) 

<Si = {a\ai — c4 a 2), (3b) 

5 2 = {a\a 2 + a 2 ai), (3c) 

5 3 = -i(a\a 2 - a\a{) ; [a i: a]] = (3d) 

The expectation values (Si) , i = 0,1,2,3 can be identified with the Stokes parameters 
Si , i = 0,1,2,3. 

As is well known, the above operators furnish the Schwinger representation || of the SU (2) 
algebra 

[Su Sj] = 2ie ijk S k , k = 1, 2, 3 (4a) 

[So A] =0, * = 1,2,3 (4b) 

S (S + 2) = S 2 X +S^+Sl (4c) 

The connection to the SU (2) algebra is very useful and many of the concepts devel- 
oped in connection with spin squeezing J| have been taken over to discuss polarisa- 
tion squeezing HJ^, p!0| , p!lJ| . The Stokes parameters Si ,i = 0,1,2,3 involve four averages 
(a\ai), (a 2 a 2 ), (a[a 2 ), (a 2 a\). Further, the variances 

V ij = ±(({S i ,S j })-{(S i ),(S j )}) ; i,j = 0,---,3,i<j, (5) 
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can be expressed in terms of the nine normally ordered field correlations 
{a\a\aiai), (a} 2 a} 2 a 2 a 2 ) , (a\a\a 2 a 2 ), (a2 a 2 a i a i)> (a\(i2 a i a 2) , (a\a[aia 2 ), (02020102)1 (<A a 2 a i a i) 1 
{a\a} 2 a 2 a 2 ) . Given the circumstance that, in an experiment, one typically measures average 
intensities and intensity variances, the question one is confronted with is whether it is pos- 
sible to deduce the relevant field correlations by measuring average intensities and intensity 
variances of suitably defined modes. The answer to this question is provided by a theorem 
due to Mukunda and Jordan [I2| in a more general context. This would be used in the 
following to show how all Vy can be essentially obtained from intensities and intensity fluc- 
tuations. We note that very recently Korolkova et al have proposed a method to measure 
the diagonal quantities Va. An implementation of this proposal has been achieved by Bowen 
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Before we proceed further, we note that for an important special class of the states of 
the radiation field, such as the squeezed coherent state |cki, Qf2> C)j ^ ne correlations can be 
derived in closed analytical form. In this case it is well known that the Wigner function as 



well as the Q — function are Gaussian | l3fl . Thus the fourth order expectation values can 



be expressed in terms of second order expectation values. If we define displaced modes dj 
<2j — CKj, then we have for example 



(did 2 44> = {did 2 ){d\d\) + {d l d\){d 2 d\) + (d 1 d\)(d 2 d\) (6) 
= (d 1 d 2 )(d\dl) + (d 1 d\)(d 2 dl), (7) 

since for the squeezed state {d\d\) = 0. It should be borne in mind that the special case of 
coherent states is obtained by setting ( = 0. Clearly if one knows, a priori, that the field is in 
the squeezed state then all the fluctuations and correlations among Stokes parameters can be 
related to second order expectation values. For arbitrary non classical states characterised 
by non-Gaussian Wigner function, the situation would be much more complex. Another 
simple case is obtained if the fluctuations of the field around the mean value are small. In 
such a case Sk — (Sk) can be approximated by terms linear in dj and df. Thus correlations 
Vij are expressible in terms of the expectation values of quadratic forms. The latter, in many 
cases reduces to variances in quadratures 

As noted earlier, The measurement scheme proposed here is based on a work of Mukunda 
and Jordan ]12| which, In the present context, involves use of unitary transformations defined 









C:) 













by 

cos 9 e 1 ^ sin 9 \ , . 

-e^sinfl cosfl J' ^' 

and measurements of the intensities (b\bi)(i = 1,2) and intensity correlations like 
(b\b\bib2) for different settings of 9, (f>. Thus, for the mean intensity, we have 

(b\(u)bi(u)) = cos 2 9{a[ai) + sin 2 ^^^) 

+ sinO cos e(a\a 2 ) + e~ i4 sin 9 cos 9 {a\ax ). (9) 

We can, therefore, determine (a|ai), (0202), (0102), (c4 a i) from the knowledge of 
{b\(u)b\{u)) , for instance, for the four settings corresponding to (9, 0) equal to 
(0,0), (tt/2,0), (vt/4, 0), (vr/4, n/2). This is, in fact, identical to the procedure adopted 
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for measurement of classical stokes parameters |L|. 



Next we turn to the nine field correlations which arise in the expressions for the variances 
of the Stokes operators. On selecting the unitary operator in (8) as 

/ cos6> sin# \ , 1rA 
U= ^-sin6 cosO )' (10) 



we find the relation 



(b\(u)b\(u)bi(u)bi(u)) = (a\a\aiai) cos 4 ^ + (a 2 a 2 a 2 o-2) sin 4 # 

+ (4(a{a 2 aia 2 ) + (a\a\a2<X2) + (a 2 a 2 aiai)) cos 2 9 sin 2 9 
+ 2((a{a|aia 2 ) + (a\a\ai<ii)) cos 3 #sin# 

+ 2((a{a 2 a 2 a 2 ) + (ala^a^)) cos 6* sin 3 6*, (ll) 

which for five different values of 9 allows us to determine the terms with different dependence 
on 9 in terms of the measured values (b\(u)b\(u)bi(u)bi(u)) for the chosen 9 settings. Next 
taking the unitary operator as 

/ cos# ism9\ /iriN 

u = , ■ • a a ■ ( l2 ) 

\+ism9 cos 9 J 

we get the relation 

(b\(u)b\(u)bi(u)bi(u)) = (a\a\aiai) cos 4 9 + (a2 a 2 a 2«2) sin 4 ^ 

+ (4(a|a2aia2) — (a|a|a 2 a2) — (a^a^aidi)) cos 2 9 sin 2 9 

+ 2i((a\a\aia2) — (a\a\aiai)) cos 3 #sin# 

+ 2z((a|a2a 2 a2) — ( a 2 fl2 a i a 2)) cos # sin 3 9. (13) 

Here, three suitably chosen values of 9 are sufficient to determine the terms with different 
dependence on 9 since the first two terms were determined previously. Finally, choosing 

U = T2\-e~^ 1 J' (14) 

we obtain 

(b\(u)bl(u)bx(u)b 2 {u)) = -(ajojaiai) + -(44a 2 a 2 ) 

— -((a[a[a2a2) — (alalaidi)). (15) 

Thus the nine choices of u above permit us to determine the nine field correlations listed 
after Eq.(5) and hence all the variances in the Stokes operators. Besides the variances, the 
same nine correlations determine the correlations among Stokes operators. For example 
the last term in Eq.(15) represents the normally ordered correlation ( ° 52S3 ° ) . Further 
note that (11) [(13)] for 9 = it/A and 9 = 37r/4 when added and subtracted would yield 



, o SqSo o ) + ( o S2S2 o ) 



( o 'S'o'S'o ° )-( ° S2S2 o ) 



and (:s Q s 2 :) [(:s s 3 i)). 
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The unitary transformations u involved in the measurements above are all of the form 



u(6,4>) = -id, - a a > ( 16 ) 

v ' \ — e l<p smv cos 9 J 

with taking values 0, 7r/4, tt/2. These transformations can most conveniently be realized 
using the universal £77(2) gadget |T5| , p!6| which consists of two quarter wave plates and one 



half wave plate mounted coaxially. Taking the axis of the gadget to be along the z axis and 
the slow axes of the quarter and half wave plates in the x direction, the settings for the 
SU{2) gadget in the Q-Q-H configuration which realize the required transformations are 

u(9,0) = Qtt/4 Qtt/4 -ff-Tr/4+0/2, (17a) 
u(0,7r/2) = Q n / 2 Q(e+7T/2) H e /2, (17b) 
w(7r/4,7r/4) = Q ff/4+(tail -i v2)/a Q^n/n+it^ V2)/2 H */i2, (17c) 

where the the subscripts to Q and H give the angles by which the quarter and the half wave 
plates have to be rotated around the z axis. The matrices for the quarter wave and half 
wave plates are given by 

. / cos 20 sin 2(p \ , . 

^ = Z Uin20 -cos 20 J' (18a) 

^=4f( cos2 t; z ^1 ■)■ ( 18b ) 

v V sin 20 —cos 20 — ij 

After 'rotating' the a± and 02 into b\{u) and 62 (u) modes with the help of the SU(2) gadget 
one can separate the bi(u), bz(u) modes with the help of a polarising beam splitter and 
carry out the required intensity and intensity-intensity correlation measurements to obtain 
complete information about the Stokes parameters and fluctuations thereof. 

To conclude, we have shown how a combination of quarter wave and half wave plates can 
be used to measure not only the Stokes parameters but also all the quantum correlations 
and fluctuations in the Stokes parameters. 
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